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ON THE FUNCTIONAL EQUATION aU + C*( X U) = F 

PHILIPPE RYCKELYNCK* LAURENT SMOCH* 

Abstract. We study in this paper the functional equation 

cm(£) + C * (xu)(t) = f (0 

where a E C dxd , u, f : R — >• C d , u being unknown. The term C * (xu)(t) denotes the discrete 
convolution of an almost zero matricial mapping C with discrete support together with the product 
of u and the characteristic function x of a fixed segment. This equation combines some aspects 
of recurrence equations and/or delayed functional equations, so that we may construct a matricial 
based framework to solve it. We investigate existence, unicity and determination of the solution to 
this equation. In order to do this, we use some new results about linear independency of monomial 
words in matrix algebras. 

Key words. Difference equations, generalized derivatives, matricial polynomial functions, ma¬ 
trix sequences 
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1. Introduction. Let [fo,t/] be some interval of time and x its characteristic 
function, and e > 0 a fixed time delay. Let ( Ck)kez be an almost zero sequence of 
matrices in C dxd where d denotes the “physical” dimension. We define next C : R —»• 
, , , , I c_ t / e if t £ eL 

C dxd by C(t) = < . We denote by Id the identity matrix of size d 

I 0 otherwise 

and we fix throughout a matrix a £ C rfx<i . We have the discrete convolution product 
over eL\ 

(u*v)(f) = ^ u(s)v(f — s), t £ R 

where u and v are two mappings defined in ]R with matrix values of format (dxd,dxl) 
or (1 x d, d x 1). 


This paper is devoted to the study of the following functional equation 

cm + C*(\u) = f (1.1) 

which may be rewritten equivalently as 

au(t) + ^2 Cfcu(t + ks)x{t + ke) = f (t), (1.2) 

fcE Z 

where u, f : R —> C d , u being unknown. The sum in the left hand-side of m will 
be denoted throughout the paper by Du = C*(yu), where the operator □ may occur 
in the context of discrete calculus of variations (see p4j [5] ). Equation (11. ip may be 
thought as a discrete version of the differential equation au(t) + u'(t) = f(t) where 
u' has been replaced by □. An analogous functional equation appears in the field of 
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time-scale calculus and has led to many works, see for instance mm and the literature 
therein. However, we do not require here that the operator □ behaves as a derivative. 


We define TV as the largest integer such that cat ^ 0 or c-n ^ 0. For example, 
when TV = 1, equation m rewrites as 


0 


f(f) — au (t) = < 


ciu(t + e) 

ciu(f + e) + c 0 u(T) 

ciu(f + e) + c 0 u(f) + c_iu(t 

c 0 u(t) + c_iu(t - e) 

c_iu(f — e) 

0 


e) 


if t < to — e 
if To — £ < T < To 
if t 0 < t < t 0 + e 
if to + e < t < tf — £ 
if tf — e < t <tf 
if tf < t < tf + e 
if tf + £ < t 


(1.3) 


In general, the equations ed may be thought as a mixture between recurrence 
equations and delayed functional equations. 

We focus in this paper on the following two problems. What are the analytical 
properties of functions of the shape Du and is the functional equation ED well-posed? 

The paper is organized as follows. In Section 2 we provide some analytical prop¬ 
erties of the functions Du. Section 3 presents some technical lemmas for matricial 
linear recurrences. Section 4 is concerned with the non-vanishing property for en¬ 
tries and determinants of matricial polynomial functions. It ensures that the linear 
compatibility equations extracted when solving ED are Cramer for the specific case 
TV = 1. We prove in Sections 5 and 6 that ED is well-posed if a is invertible, with 
proper accuracies, according to the respective cases TV = 1 and arbitrary TV. The 
case a = 0 is also investigated, and stands for the characterization of the range of the 
operator □. At last, Section 7 gives some concluding remarks and perspectives. 


2. Analytical properties of the functions Du. In order to comprehend equa¬ 
tion ED and present some robust methods for solving it, we provide first some fea¬ 
tures of the operator □. 

Theorem 2.1. Let i^eNU {oo}. 

1. The operator □ maps the space of the functions u : R —> C d of piecewise- 
C 1 ' regularity, to the space of compactly supported functions v : M —► <C d of 
piecewise-C 1 ' regularity, such that Supply) C [To ~ TVe,T/ + TVe]. 

2. If u : K —> C d has p points of discontinuity on R, then Du has at most 
(4TV + 2) + (2TV + 1 )p points of discontinuity on R. 

3. If cjv or C-n is invertible then the kernel of □ consists of the functions which 
vanish on [fo,t/]. 

4- If Cn or C-n is invertible then, for all u : R —> C , Du is measurable, 
(respectively integrable, of piecewise-C 1 ' regularity) on R if and only if u is 
measurable (respectively integrable, of piecewise-C 1 ' regularity) on [Tq,T/]. 
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Proof. The main ingredient of the proof is the explicit formula for Du(t) provided 

in 0: 

0 if t < to — Ne 

Civu(f + Ne) if to — Ne < t < to — Ne + e 

CAr_iu(f + Ne — e ) + cjvu(f + Ne) if to — Ne + e < t < to — Ne + 2 e 


cjv_ p u(f + Ne — pe) + ... + cjvu(t + Ne) 


< c_jvu(f — Ne) + ... + cjvu(f + Ne) 


if p S (0,..., 2 N — 1} and 

to + (p ~ N)e <f<to + (p + l — N)e 


if to + Ns < t < tf — Ne 


( 2 . 1 ) 


C-N\i(t — Ne) + ... + c p _atu(£ — Ne + pe) 


if p G {2 N — 1,..., 0} and 
tf + (N — p — l)e < t < tf -b (N — p)e 


c_Aru(t — Ne) + ci-Nu(t — Ne + e) if tf + Ne — 2 e < t < tf + Ne — e 

C—Nu(t — Ne) if tf + Ne — e < t < tf + iVe 

0 if tf + TVe < t 

Assertion 1. is an obvious consequence of formula (12.11) . Next, let u : R —> C d , then 
formula (12.11) shows that the (AN +2) points to + ke and tf + ke, \k\ < N may be points 
of discontinuity of Du. Moreover, if tg is any point of discontinuity of u, the values 
tg + fee, |fc| < A r , give rise in general to points of discontinuity of Du. Lastly, since 
there does not exist any other point of discontinuity of Du, assertion 2. is proved. 
Let us prove assertion 3., i.e. the conditional injectivity of □. We first notice help to 
m that if u vanishes in [to,t/], then all the terms occuring in (12.11) are equal to 0 
so that Du is zero everywhere. Conversely, suppose that det(cAr) ^0. If u G ker(D), 
second row of formula m shows that u = 0 in [t 0 ,to + £[- Next, by using this result, 
the third row of (EH) shows that u = 0 in [to + e, to + 2e[ and so on. Then we may 
prove easily that u = 0 in [t 0 , to +e[U[to + £, to + 2e[U ... U]t/ — e, tf]. We may proceed 
similarly when det(c_Ar) 0 by starting from the last but one row of (12.11) . As a 
consequence, if u, v : R —> C d , we get 

□ u = Dv u = v in [to, t/] 

which proves the third assertion. We proceed in the very same way to prove that 
measurability or integrability of u on [io,t/] is equivalent to the same property for 
□ u on R. □ 


3. Preliminary lemmas for matricial linear recurrences. As we shall see in 
Sections 5 and 6, solving EH) queries to deal with matricial nonstationary recurrences 
of the shape 


w„+i = M„w„ + g„, n G N, 


(3.1) 


where M„ is a s x s matrix and w„,g„ are vector sequences in C s , s being a fixed 
integer. Since non-commutative products occur in the explicit formula, let us denote 


3 




by 


]jM fc =M m M ro _ 1 ...M 0 

k =0 

this left-side product. 

Lemma 3.1. The solution to the recurrence mu is given by 


w „ 




g k- 


(3.2) 


Proof. The cases n = 0 and n = 1 are true and use the conventions of the 
empty sum and empty (left-)product. The proof for n > 1 is straightforward by using 
induction. □ The particular case when M is independent on n is of some importance 
and gives rise to the solution to stationary matricial recurrence of the shape m as 
follows: 

n— 1 

w n = M"w 0 + ^M ra - fe - 1 g fe . (3.3) 

fc=o 

In the following two lemmas, we deal with sequences of non-commutative matricial 
polynomials. 

Lemma 3.2. Let d > 2, j3,j € A4d,d(C) and ( 5 ' n ) and (<5") be the sequences of 
matrices defined by the recurrences 

S'n+l = P S h + <5" + l = O + <^'-l7> n > 1 (3.4) 

with 5' 0 = Sq = 0 and = 5" = Id- Then, for all n>2, we have 

S' n = <5" = ... (3.5) 

where the sum is extended over all the multiplets (mi, m 2 ,..., m n ) € N" with 

[5] 

irij >1 if 2 < j < n — 1 and J^(m 2 i-i + 2m2i) = n — 1. (3.6) 

2=1 

The number of monomials in the formula 1^.51) is the Fibonacci number T n . 

Proof. First, let us mention that formula (13.51) is true for n = 0 and n = 1 by 
using the empty sum convention and the initial conditions. Let us prove now the 
formula (13.51) for the sequence (5' n ). Suppose that the result (13.51) holds for n and 
n — 1 and that 5' n _ 1 and 5' n are of lengths T n -1 and T n respectively. Then, help to 
(13.51) and (13.41) we see that d' n+l is the sum of the words excerpted from the sum 5' n 
multiplied on the left by [5 and those from the sum multiplied on the left by 7 . 
So we have proved by induction that S' n has the shape 

* = £ e ( n, mi , m 2 ,...)/9 mi 7 m 2 /? m 3 7 m4 ■ • ■ 
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where the coefficients e^ n , r n 1 , r n 2 ....) are convenient positive integers independent of /3 
and 7 . Let us note that if /? = 7 = Id then, obviously, S' n = T n x Id so, the number 
of monomial words in the sum (13.51) is T n . Let us prove now that all the coefficients 
are in. fact equal to 1. Indeed, we see that all the T n ~ 1 words coming 
from 8' n _ 1 are distinct from the J~ n words arising from 8' n since the first ones begin 
with /3 while the second ones start with 7 . All the words S' n+1 occuring in m being 
distinct, their cardinality is T n ~\ + T n = T n +i- So far, we have proved the result 
(ESI) for S' n and the proof for 5" is the same. □ In the following, we set 8 n = S' n = <5", 
Vn G N. 

Lemma 3.3. Let d> 2, (3,^ G Afd,d(C), (gG (C d ) N , and let (w n ) n G (C d ) G 
N the vector sequence satisfying for all n > 1 the following two-order matricial recur¬ 
rence equation: 

w„+i = /3w„ + 7 W„_i + g n . (3.7) 


Then for all n € N* we have 

n —1 

w n = (5 n Wi + ^n-l7 w 0 + ^2 tin-kgk- (3.8) 

k=1 


Proof Let us define M = 


0 7 

Id 0 


We may easily prove by induction that 


M" = 


8n +1 

8n dn—l'y 

(S n ). Next, (13.71) may be rewritten as 

'w „ +1 


, for all n G N*, by using the recurrence (13.41) for the sequence 


By setting = 


w„ 

w„_i 


Sn = 


= M 

Sn 

0 


w„ 

W„_l 


and by using a slightly adapted version of 

formula (13.3D . we get w ra = M n_1 wi + J2k=i M n ~ x ~ k gk■ Now, the formula (13.81) is 
obtained by reading the first row of the previous formula. □ 


4. Linear independency of words in matrix algebras. We shall obtain 
a result which ensures that the foregoing linear compatiblity equations which arise 
when solving the functional equation EH) are Cramer for almost all choice of matrices 
a, c_i, Co, Ci G C dxd if N = 1. 

Let us call a matricial polynomial function of two matricial indeterminates (£, £) 
an expression of the following shape 


<hc, o = E ^ lim2 ,...c mi r 2 c ro3 r 4 ■ ■ • (4.1) 

with coefficients in C. As an example, we may cite the sequence of matrices occuring 
in Lemma 13.21 Given such a function, for all (p, q ) G N 2 , we denote by K pq the sum 
of all the coefficients K mi ,m 2 ,... such that mi + m 3 + ... = p and m 2 + mi + ... = q. 
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Similarly, for all (p, q,r, s) £ N 4 , we introduce the sum K p q r s of all the coefficients 
mi x m 2 x ... x Kmi,m 2 ,... such that mi +TO 3 +... = p, TO 2 +TO4 +... = q, and where r 
and s count the occurences of ( and f respectively in the monomial of HD associated 
to the multiplets (mi, m 2 ,...). We must have necessarily \r — s| < 1. We shall say 
that the function hd is generic when at least one sum K pq and one sum K p ^ r , s 
are nonzero. For instance, the matrices S n occuring in (13.51) are the specializations at 
(/?, 7 ) of generic polynomials, since the coefficients are either 1 or 0 . 

Let us give now the main result of this section. 

Theorem 4.1. Suppose we are given l generic matricial polynomial functions 
an d us se t V’i k j (C) £) = £))ij- Then, the set of couples (C>£) 

such that 

\/i,j,k, V’|j ) (C>£) 7 ^ 0 , with 1 < i,j < d, 1 <k <1, and 
Vfc, det(£,£)) 7 ^ 0, with 1 < k < l 
consists in an open dense subset of (C dxd ) 2 . 

Proof In order to prove this result, we first show that if ip(C,f) is the generic 
matricial polynomial function defined by (ED, then all its entries ipij (£, £) are nonzero 
polynomials w.r.t. their 2d 2 indeterminates. The proof goes as follows. Let us 

specialize if at matrices of the shape 

( x z 0 ljd _ 2 \ / y t 0 M _ 2 

0 x 0 M _ 2 J and £ = I 0 y Q14-2 

Od-2,1 0d_2,l Id-2 J \0d_2,l 0d_2,l Id-2 

for convenient x,y,z,t £ C. We easily get 

/ if(x,y) %p(x,y,z,t) 0 M _ 2 \ 

^(C,£)=l 0 ip(x,y) 0 1)d -2 ( 4 - 2 ) 

V 0d_2,l 0d_2,l V>(1, 1)^-2/ 

where ip is the function HD evaluated at complex numbers and is equal to ^ Kp,qX p y q 
while ip is the polynomial K P , q , r ,sX p ~ r y q ~ s z r t s . Due to the hypotheses, the sets 
of couples (x,y) or quadruplets ( x,y,z,t ) such that ip(x,y) 7 ^ 0 or ip(x,y, z,t) 7 ^ 0 
respectively, are open everywhere dense subsets in C 2 and in C 4 respectively. Now, 
in order to deal with the zero-like entries of (14.21) . we use the following property of ip. 
We have 

V’(S'“ 1 CS5S' _1 £S') = 

where g is any invertible matrix in C dxd _ Especially, if we choose g as a permutation 
matrix p a , we obtain 

The symmetric group e d acts transitively on the set of couples (i,i) with 1 < i < d 
and transitively on the set of couples (i,j) with 1 < i, j < d and i 7 ^ j. Thus, since 
ipip and ipi t 2 are nonzero polynomials, we may claim by using suitable specializations 
of (C,£) that every other entry ipij of ip{<p,0 gives rise to a nonzero polynomial. 
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Let us prove now that if ip is generic, then det(ip(£u ,... ,£,dd)) is a nonzero poly¬ 
nomial. Indeed, if we particularize £ and £ to be two multiples of the identity matrix, 
i.e. (Cj£) = (pldicld) with p and cr in C, we get the polynomial 

det {ip{pld, oh)) = ('^2Kp,qP P o q \ = V>(P> a ) d 

\ p,g J 

which is nonzero in C[p, cr] since there exists at leat one coefficient K pq 7 / 0. 

Collecting all these informations, we may work with open dense sets in C 2d 
instead of dealing with nonzero polynomials. Let us consider for all k £ {1,..., 1} and 
i. j £ { 1 ,..., ci} the closed algebraic varieties of codimension 1 

= {(C.O e (c dxd ) 2 /^?(Cn, • ■ • ,Zdd) = 0} 

and 

■TdS ={(C,0 G (C dxd ) 2 /det(^ (fe) (Cn,--.,^d)) =0}. 

Since each hypersurface is meager, the union (U ijkJ'-j ' 1 ) [J(UfcJ r ^) is a closed meager 
subset of (C dxd ) 2 . Thus, its complementary is an everywhere open dense subset of 
(C dxd ) 2 . This ends the proof of the result. □ 

Remark 4.1. As an example, for all integer m, the set of couples (/3,7) £ (C dxd ) 2 
such that 


det(/3) 7 / 0, det( 7 ) 7 / 0, /?7 7 ^ 7/3 and det( 6 n ) 7 / 0, Vn £ {0,..., m}, (4.3) 


is an open dense subset of (C dxd ) 2 . In the very special case when the matrices /3 ,7 
are multiples of Id, i.e. /3 = /3/d and 7 = 7 Id with /3 and 7 in C, we easily find help 
to (13.41) that 

( P+y/ 18 2 +47 


Sn = 


Id- 


P—y/p 2 +i : i 

2 J V, 

The set of couples 0,j) such that is not °f the shape exp (^yp 1 ) is open 

dense in C 2 and, for those couples, we have det(<5„) 7 / 0. 

Remark 4.2. The set of generic polynomials is an open dense subset in the vector 
space of matricial polynomial functions of given degree. 

Remark 4.3. Let us mention how to proceed if we deal with non-generic polyno¬ 
mials. We may use in this case the so-called defect Theorem of P.M. Cohn, see |3J 
Theorem 9.6.1, p. 283]. If /C is a commutative field and H is an alphabet of non com¬ 
muting variables, we may denote by IC((Cl)) the algebra of non commutative formal 
series in these variables. Then the defect theorem states that if are elements 

of K.((H)) without constant term, satisfying a non-trivial relation <^ 2 ) = 0 for 

some non commutative series 4> in two variables then <pi , ip 2 commute. This theorem 
allows to exhibit at least one nonzero polynomial entry to ip(C,,£,), but the conclusion 
is not as so accurate than in Theorem 14.II 
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5. Solving explicitly the functional equation when N = 1. Let us show 
now how to solve the functional equation m when N = 1 by using the preceding 
results. We restrict ourselves to the quite opposite cases a invertible or a = 0. Let 
us define M > 1 as the integer part of tf ° . 

Theorem 5.1. If N = 1, for all quadruplets (c-i,co,ci,a) in an open dense 
subset of (C dxd ) 4 , the equation 

cm + C * (yu) = f 
admits one and only one solution. 

Proof. We may assume by density that c\ is invertible. If M = 1, the result 
follows easily from the inspection of (11.31) . 

The main part of the proof consists in solving m and is valid for all a. Outside 
[to — e,tf + e], one has f(t) = cm(f) which allows to determine u if a is invertible. 
Next, we focus on the interval [to + e, f/ — e]. By using (11.21) . we get the recurrence 
equation: 

u(r + (n + l)e) = + c 0 )u(r + ne) - c_iu(r + {n — l)e) + f(r + ns)), (5.1) 

Vr 6 [to, t 0 + s[ and n > 1. So we see that (15.11) is a recurrence of the shape (13.71) . 
Using Lemma 13.31 with 


/3 =-c 1 1 (a +c 0 ), 7 =-c 1 1 c_i, (5.2) 

and w n = u(r + ns), g n = cf 1 f(r + ns) for r fixed, we get 

n—1 

u(r + ne) = S n u(r + e) + J n _i 7 u(r) + ^ 5 n _fcC^ 1 f(r + ke), (5.3) 

fc=i 

for all instant r + ns < tf. Since r £ [to,to + e[> we have in any case n < M — 1. 

More explicitly, we shall use (15.31) in the following ranges of indices: 

t £ [t-o,tf — Me[=> max{n/r + ns < tf} = M , (5.4) 

t £ [tf — Ms, to + e[=> max{n/r + ns < tf} = M — 1. (5.5) 

We find it convenient to denote by I = £(r) the maximal value of n for which we 
may use (15.31) . given by do~H) or (15.51) . Although the formula (15.31) is true for n = 1, 
this result is tautological in this case. For n > 2, (15.31) expresses the solution to the 
functional equation m in the interval [to + 2e, t /], as a linear combination of the 
restrictions u(r) and u(r + s) of u to the two intervals [to, to + s[ and [to + s, to + 2e[. 
The coefficients of this linear relationship are independent on time and constitute the 
sequence of matrices (S„) which depend only on /3 and 7 . 

Now, let us determine the two additional unknown functions u(r) and u(r + e) by 
solving near the boundaries to and tf. In order to do this, let us write explicitly 



the four remaining equations: 


if t G [to — £,to[, f(£) = cm(i) + ciu(f + e), (5.6) 

if t G [to, to + e[, f(t) = (a + c 0 )u(t) + Ciu(t + e), (5.7) 

if t G]t/ - e,tf\, f(t) = (a + co)u(t) + c_iu(t - e), (5.8) 

if t£]tf,tf+e\, f(t) = au(t) + c_iu(t — e). (5.9) 


We note that the restriction of u to [to — e,to[ is undetermined. Let us call this 
function tp. Next we solve m and as follows, by adapting slightly the range of 
the time variable. If t G [to, to +e[, (15.611 provides 

u(t) = cj" 1 f(t — e) — Ci 1 aip(t — e). (5.10) 

If t G [to + e, to + 2e[, (15.71) and (15.101) give 

u(t) = c(" 1 f(t — e) — cj” 1 (a + co)c(" 1 f(t — 2e) + c)" 1 (a + co)ci l a<p(t — 2e). (5.11) 

The formulas (15.101) and (15.111) give the two additional restrictions of u needed in the 
expansion <53*. We note that these restrictions, at the time being, depend on the 
auxiliary restriction tp of u. 

Now it remains to solve (15.81) and (15.91) . First, the two functions occuring in the 
r.h.s. of the equation (15.81) are already known and may be expressed through (15.31) . 
The equation (15.81) may be thought as a constraint either on ip or on the l.lr.s. f 
of (11.11) . We may observe that tp(r — e) occurs two times in u(£) and two times in 
u(£ — e). So, the coefficient we are looking for is equal to 

9 = ((a + cq)8i + c_i<^_i)coeff of p in u(r + e)] 

+((a + co^-i + c_i<^_ 2 ) 7 [coeff of <p in u(r)], 
where l = M or i = M — 1, depending on the location of r w.r.t. tf — Me. Direct 
inspection of these coefficients from (15.101) and (15.111) yields the explicit expression of 
the coefficient of <p in (15.81) : 

9 = ({a + c 0 )S e + c_i^_i)(c)( 1 (a + c 0 )cj" 1 Q!) 

+((a + c 0 )5e-i + c_id f _ 2 )7(-cr 1 a) (5-12) 

= [((a + c 0 )5e + c_i^_i)c)( 1 (a + c 0 ) - ((a + c 0 )^_i + c_i^_ 2 ) 7 ]c[" 1 q:. 

Now, in order to simplify the matrix 9 , we use the formulas <53* to eliminate a + cq 
and c_i and we get 

9 = ci((/38i + 7<^_i)/3 + (fiSi-i + 7^_ 2 )7)c]" 1 a. 

By using three times the recurrence (13.41) . we obtain 

9 = cid^+ 2 c 1 1 a. 

At the time being, no assumption has been made about a. From now on, we consider 
the case when a is invertible and chosen in such a way that det(i5f+ 2 ) ^ 0 (see (15.21) 
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and Remark H~T1) . Therefore, det(d) ^ 0. Hence, (15.81) may be solved w.r.t. ip(r — e). 
As a by-product of the previous results, <p(r — e) is a linear combination of the values 
of f on the set {r + ks\k\ < N}. It remains to solve equation (15.91) on ]t/,t/ + e]. 
We get simply u(t) = a _1 (f(t) — c_iu(i — e)) for all t g]t/,t/ + e], where the r.h.s. 
has already been determined. Then we have proved that, for a generic quadruplet 
(c_i, Co, Ci, a) and for all f : R — > C d , there exists one and only one solution to the 
functional equation m- □ Next, let us consider the case when a is zero, which 
stands for the characterization of Im(D). 

Theorem 5.2. For all triplets (c_i,co,ci) in an open dense subset of (C dxd ) 3 
and for all g : [fo,f/] —> C d , there exists one and only one extension f of g to K. such 
that Supp(i) C [to — s,tf + e] and the functional equation f = Du admits one and 
only one solution in [to,t/]. 

Proof. We may assume by density that c\ and c_i are invertible. Let be given a 
function g : [t 0 , f/] —> C d and let us construct an appropriate extension f : K. — > C d of 
g as follows. We set 


0 

if t < to — £ 

p(*) 

if to — e < t < to 

s(t) 

if to < t < tf 

q(*) 

if tf < t < tf + e 

0 

if t > tf + e 


for some auxiliary functions p and q with values in C d . The proof consists in showing 

that f lies in the range of □ if and only if p and q are uniquely determined by g. We 

keep the notations of the previous proof. When a = 0, the formulas (15.11) to (15.111) 

still hold, we have 9 = 0 and the auxiliary function <p does not occur anymore. Let 

us solve the two equations (15.81) and (15.91) w.r.t. the two restrictions p and q of f. To 

write these equations, we introduce the operators 

t -1 

A*f(r) = ^2 f(T + ke). 

k =1 

Let us notice that these operators may be thought as discrete integral operators. 
Plugging formula (15.31) . with the maximal value l allowed for n, in the equation (15.81) . 
we get for all t G]t/ — e,£/] 

f (*) = co(&u(t + £ ) + <^-i 7 u(t) + A^f(r))+ 

c_i(5<_iu(r + e) + <^_ 2 7u(t) + A^_if(r)), 
t and r being connected as previously by the requirements that r £ [to, to + £ [ an d 

g N. Remembering the values of /3 ,7 given by (15.21) and the recurrence (13.41) . we 

get: 

f(t) = -ci(^+iu(r + e) + ^qu(r)) + c 0 A f f(r) + c_iA*_if(r). (5.13) 

Since we may respectively rewrite (15. 101) and (15.lip as 

u(r) = cf 1 p(r) and u(r + e) = cj" 1 f(r) - cf 1 c 0 cf 1 p{T), 
we see that the coefficient of p = f(r — e) in (15.131) is given by: 
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(<^+i c i lc o ~ ^7) c i 1 = -^+2C X x . 

We assume the triplet (c_i, c 0 , Ci) lies in the open dense subset where det(^ +2 ) 7 ^ 0 
and thus, we have shown that the restriction p exists and is unique. This being done, 
the equation (15.81) yields 

f (t + e) = q(i) = c_i(d £ u(r + e) + <^_i 7 u(t) + A^f(r)), 

t and r being connected as before, which is well-determined. As a consequence, the 
formula (15.31) gives rise to one unique function u on [to, t/] such that Du = f and this 
ends the proof.D 


6 . A matricial-based framework when N is arbitrary. Let us use a matricial- 
based framework for solving the functional equation for arbitrary N. In this section 
we assume throughout that cjv and a are invertible, mainly for the ease of exposition. 
We shall say that the (2A + 2)-tuple (c_at, • • ■, Co,..., Cat, a) is generic in (C dxd ) 2Ar+2 
provided the foregoing system (16.111) is Cramer. This system (16.111) consists in N equa¬ 
tions w.r.t. N unknowns, generalizing the auxiliary function (p occuring in Section 5. 
Its coefficients are matricial polynomial functions w.r.t. cf, 1 (cq + a) and the 2 N — 1 
matricial indeterminates cffCk, k g {— N ,..., N — 1} and k 7 ^ 0, which in a sense 
generalize m- 

Theorem 6.1. If the (2N + 2)-tuple (c-N, ■ ■ ■ ,co,... ,Cn,cx) is generic, the func¬ 
tional equation m admits one and only one solution. Similarly, we shall say that 
the (2 N + l)-tuple (c_at, ..., Co,..., cn) is generic in (C dxd ) 2Ar+1 provided the fore¬ 
going matrix 0 defined by (16.121) is invertible. This matrix is associated to a linear 
system of N equations w.r.t. N unknowns pi,... pjv, which are the restrictions of f 
to parts of [ffi — Ale, to[ in order that f lies in the range of the operator □. 

Theorem 6.2. If the (2N + 1)-tuple (c_jv, .. ., co,..., cjv) is generic, for all 
g : [to,t/] —> C d there exists one and only one extension f : R —» C d such that 
Supp( f) C [to — Ne ,t f + Ne\ and the functional equation f = Du admits one and only 
one solution in [to,t/]. 

Proof. (Of both theorems.) 

Let u : R —> C d be a solution to equation m, related to A, B and U defined over 
R as follows: 



( c N f (t) 

0 


( u{t + (N - 1 )e)\ 


B(i) = 

e c 2dN , u(t) = 

u(f) 

g c 2dN 


V 0 ) 





V u(t - Ne) J 


( 6 . 1 ) 
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A(i) = 


/-c Ar 1 c M _ix(t + (Af-l)e) ... -e^ 1 (a + cox(t)) -c^ c^ N x{t ~ Nc)\ 

I d 0 ... ... 0 


0 I d 0 J 


■ (6-2) 


The matrix A(t) G C 2dNx2dN j g i oca j]y constant w.r.t t. Especially, when t £ 
[to + Ne,tf — Ae], all the characteristic functions of the entries of A (t) are equal to 
1 and we shall agree to denote by A the constant value of A(t) in this segment. Let 
us mention that the matrix A depends only on the 2A + 2 matrices Ci and a and is 
nothing but the block companion matrix of the functional equation (ED- 

When t< to — Ne or t>tf + Ne, we have already seen that the solution to the 
equation ED is given by u(t) = a 1 f (t). Now, if t £ [to — Ne,tf — Ae], (11.21) may 
be written as 


u(t + Ae) = c Ar 1 [f(t) - (cjv-ix(t + (A - l)e)u(t + (A - l)e) + ... 

+{a + c 0 x(t))u(t) + ... + C- N x(t - Ne)u(t - Ne))]. (6.3) 

When t £ [t 0 — Ne, t 0 [, this equation expresses u(t + Ne) explicitly as a combination 
of the A auxiliary functions = u(r — ke), 1 < k < A, which are at this 

level unknown. Obviously, the functions 1 < k < A, stand for the last A 

components of U(r). By using definitions (16.11) and (I6.2L recurrence (16.31) may be 
rewritten as 

U(t + e) = A(t)U(t) + B(t), Vt G [t 0 - Ne, tf - Ae]. 

The preceding recurrence and Lemma 13.II yield the following formula 


U(t + ne) = C„(t)U(t) + F n (t) (6.4) 

for to — Ne < t < t + (n — l)e <tf — Ne, where 


C„(t) = 


A(t + ke) if t + (n— l)e < to + Ne, 

k=0 

A n ~ N ( A(t + ke) | otherwise 


, fc=o 


and 


n —1 / n —1 \ 

= E II A (t + te) B (t + ke). 

fc=0 \^=fc+l / 


(6.5) 


( 6 . 6 ) 


Focusing on the first row of (16.41) . i.e. the first component u(t + (n + A — l)e) of 
U(t + ne), we emphasize on the fact that the unknown function u(t) may be expressed 
as a linear combination of the delayed functions i(t+ke), k £ Z through F„(t), as well 
as the auxiliary functions </ 3 fe('r), 1 < fc < A, through U(t). So, as a first conclusion 
of these calculations, u(t) is well-determined in the whole interval ] — oo,t/], as a 
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combination of the auxiliary functions iflk- 

When t € [tf — Ne,tf + IVe], we may deduce from (11.11) and (12.11) the two systems of 
N functional equations which are analogous to (15.81) and (15.91) . i.e. 

' f(4) = au(t) + c_jvu(i — Ne) + ... + cjv-iu( 4 + (JV — l)e) in ]tf — Ne, tf — (JV — l)e] 
f(4) = a-u(4) + c_atu(4 — Ne) + ... + cjv_ ju(t + (JV — 2)e) in ]tf — (JV — 1 )e,tf — (JV — 2)e] 

f(4) = au(4) + c_atu(4 — Ne) + ... + eou(t) in ]ty — e,tf] 

and 

' f (4) = au(4) + c_jvu( 4 — JVe) + ... + c_iu(4 — e) in ]tf,tf+e] 
f (4) = au(t) + c_jvu(t — Ne) + ... + c_ 2 u(i — 2e) in ]tf + e, tf + 2e] 

f (t) = au(t) + C-jyu(t — Ne) in ]t/ + (JV — l)e, 4/ + Ne] 


(6.7) 


(6.8) 


Let us consider the case when a is invertible. We convert (HD into a linear system 
for the N auxiliary unknown functions ipk- To do this, we note first that (I6.4|l . with 
n = M and t = r — Ne, writes as 


U(t + (M - N)e) = C m (t - Ne)U(r - Ne) + F M (r - Ne). (6.9) 


The first N components of U(r - Ne) which occur in the right-hand side of the 
previous equality are the auxiliary functions while the last N components are values 
of the function a -1 f at instants r— (. N+k)e , 1 < k < N. Next, by shifting adequately 
t in the system of equations HD , we get for all t G] t / — Ne ,tf — (N—l)e] the following 
equality 


( f W \ 

f(t + e) 

\f(t + (N — l)e) / 


= D 


/ u(t + (N — l)e)' 
\ u (t — Ne) 


( 6 . 10 ) 


where 


/ cn-i . ci (a + co) c_i ... c-n\ 

CN-2 ■ ■ ■ Cl (a + Co) C _1 . . . C-N 0 


D = 


Cl 


g (^Ndx (2Nd) 


\ (a + co) c_i ... c-n 0 . 0 / _ 

is a rectangular block Hankel matrix. By setting t = r + (M — N)e in (16.101) and 

using HD , we easily get 


BC m {t - Ne)V(r - Ne) 


/ f(r + (M — N)e) 
V f (r + (Af - l)e) 


DFm(t — Ne). 


( 6 . 11 ) 


By partitioning the (Nd) x (2Nd) matrix DCm(t — Ne) as (0i © 2 ), where ©i, ©2 £ 
i^NdxNd.^ we gee ^ ia j- previous linear system consists in N vectorial equations in 
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C d depending on the N auxiliary functions ip/, and is Cramer provided that ©i is 
invertible. Under the assumption of genericity, the previous condition is satisfied and 
thus, there exists one and only one solution to (16.111) . The auxiliary functions being 
determined, as a rule, the function u is well defined and unique on ] — oo,f/]. Since 
a is also invertible by assumption, we may compute the various extensions of u on 
[tf,tf + iVe] using the system (16.81) . This ends the proof of Theorem 16.II 

The proof of the second theorem proceeds in the same way. Let pj, be the unknown 
restriction of f to [to — ke, to — (k— l)e[, for all k € {1,..., N}. Similarly let us define 
as the unknown restriction of f to ]tf + {k — l)e, tf + ke], for all k £ {1, ..., N}. When 
a is zero and det(cjv) ^ 0, the system (12.11) shows that u(i) may be arbitrarily chosen 
outside [to,tf]. Indeed, the auxiliary functions ipk do not occur in the calculation of 
u. However, u is well-determined in the whole interval [to,tf] provided g is given in 
this interval and the Pfc, 1 < k < N, are already determined. In order that u satisfies 
the systems iKTTt and (16.81) . the restrictions pk and q*,, 1 < k < N, must be chosen 
adequately. Let us focus on (16.71) or, equivalently, on (16.111) . Indeed, both systems 
still hold when a = 0. However, our approach is different from previously since 
the unknowns which are the pUs are located exclusively in the term -DFm(t- Ne). 
Indeed, the above discussion shows that U(r— Ne) does not depend on the restrictions 
Pfc nor the matrix (f(r + (M — k)e))i<k< n - Let us determine the coefficient [0]^ € 
<C dxd of p j(r) in the i-th equation in the system (16.111) . To do this, we use the 
convention that for some matrix 0 £ ^(Nd)x(Nd) ^ t j ie S y m bol [0]^- denotes the (d x d)- 
block located at the (i, j)-entry of 0. Then, help to formulas (16.11) and (16.61) . we see 
that 


[®] ij 


M—l 

D A(r + (k - N)e) 

k=N—j+1 


-'N 


( 6 . 12 ) 


We assume that the multiplet (c_at, ..., Cn ) is generic, which amounts to say that 
0 is invertible. Therefore, (16.11 1) is Cramer and admits one and only one solution 

(pi,_ Pjv). Next, by using the same transformation which led to (16.101) . the system 

m may be rewritten as 


m 


\ f(t + (N — l)e ), 


( nW 


Vqjv(t)/ 


( C -1 C-2 ... C-N\ 

C—2 • • • C-N 0 


/ u(t-e) 


\ u(t — Ne ), 


V C—N 0 ... 0 / 

for all t + £■]. Since the matrix (u(t — ke))i<k<N is known, we obtain as a 

conclusion that the restrictions (pfc,qjfe), 1 < k < N, are uniquely determined by g, 
which ends the proof of Theorem 16.21 □ 


7. Final comments. Before concluding this paper, it is natural to ask if there 
exist other ways to handle the functional equation au + C* (\u) = f. Since (11.11) is a 
convolution equation, one might want to use the discrete Laplace (or Z-) transform to 
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solve it. Indeed, if the Z-transform is defined by £u(p) = exp(— pt)u(t) then 

one has the formula £(u*v)(p) = (£u)(£v) for all functions u and v on the real line, 
with compact support and values in suitable algebras. So, when a = 0, the image of 
(11.111 is £(C)£(xu) = £(/) and, thus, if £(£) is invertible, we get 

£(xu )(P) = (Efcez Ck exp(pfce)) - (Efcez exp(-pfce)f(fce)). 
Unfortunately, the determination of the original u has been transferred to a difficult 
problem of moments. Moreover, this approach does not apply to the general case 
a^O. 

Another way to proceed would be to apply the Picard’s method. Indeed, one sees 
that (ED is equivalent to 

u = a -1 f — a _1 Du = a -1 f — a _1 D(a _1 f — a _1 nu) = ... 

= a~H + EE(- 1 ) fca_fc_lnfcf - 

Thus if 11 ck 1 1| is small enough, one may hope that □ is a contraction of a suitable 
Banach or Frechet function space J-. The main difficulty in this approach is the 
requirement that the vector space T is complete, and it is related to the increasing 
number of points of discontinuity of m fc f, as k increases. 

By comparison, the approach developed in this paper is less complicated and more 
robust in the sense that the solution is given explicitly. Indeed, the central parts of 
the proofs of the four previous theorems highlight not only the well-posedness of ED 
but also an effective computational way to construct interval by interval the solution 

u. 

In scope of future work, we suggest the three following problems. First, the for¬ 
mulas (16.41) and (16.51) imply that, in the very special case when f is constant, u (t) may 
be expressed through exponential-monomial functions inside [to + Ne,tf — Ne]. So, 
in this case, the search for pseudo-periodic solutions of ED is linked to the difficult 
elimination problem of requiring that the spectrum Sp( A) of A is included in the 
unit circle of C. Second, we conjecture the following variant of Theorem 14.11 If the 
coefficients K mi ,m 2 ,... occuring in formula (14.11) are rational or even algebraic over Q 
and, in contrast, the entries (i, ? and £ij generate a subfield of C of transcendance 
degree equal to 2d 2 , then all the quantities ip-j (C,£) and det(ip^)((, do not van¬ 
ish. Lastly, we conjecture also that for all N > 2, the genericity assumption on the 
multiplets in Theorems 16.11 and 16.21 hold on open dense subsets of appropriate prod¬ 
ucts of matricial algebras. This conjecture seems natural in light of the two theorems 
presented in Section 5. 
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